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An injection of energy into the early Universe on a given characteristic length scale will result in 
turbulent motions of the primordial plasma. We calculate the stochastic background of gravitational 
radiation arising from a period of cosmological turbulence, using a simple model of isotropic Kol- 
mogoroff turbulence produced in a cosmological phase transition. We also derive the gravitational 
radiation generated by magnetic fields arising from a dynamo operating during the period of tur- 
^NJ \ bulence. The resulting gravitational radiation background has a maximum amplitude comparable 

. to the radiation background from the collision of bubbles in a first-order phase transition, but at a 

' lower frequency, while the radiation from the induced magnetic fields is always subdominant to that 

from the turbulence itself. We briefly discuss the detectability of such a signal. 
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I. INTRODUCTION 



Gravitational radiation is likely the only direct source of information about the Universe at very early times. 
Electromagnetic radiation has propagated freely only since the epoch of recombination at a redshift z ~ 1000; any 
, radiation produced at earlier times was quickly thermalized by Compton scattering from free electrons in the primordial 
plasma. Neutrinos probe to somewhat earlier epochs since they were in thermal equilibrium only until the Universe was 
around one second old, but detection prospects for the cosmic neutrino background are nil. In contrast, gravitational 
radiation was in thermal equilibrium only at temperatures approaching the Planck energy when the Universe had an 
■ age of around the Planck time. Furthermore, gravitational radiation, unlike electromagnetic radiation, propagates 
, virtually unimpeded throughout the entire history of the Universe. These properties make gravitational radiation a 
fS |- powerful probe of the very early Universe, in principle. The difficulty is of course the extremely small amplitude of 
Q ' the propagating metric perturbations. 

5^ ' The most cosmologically interesting gravitational radiation sources are stochastic backgrounds produced by some 
' event in the early evolution of the Universe. One widely discussed example is the background of tensor metric pertur- 
bations produced by quantum fluctuations during inflation Q. However, the amplitude of temperature fluctuations in 
^ the cosmic microwave background likely limits the amplitude of an inflationary gravitational wave background to be 
• I— I . undetectably small on scales amenable to direct detection (i.e. laboratory to solar system scales). Another possibility 
' is a significant background from the evolution of topological defects such as cosmic strings . Current measurements 
^ of the microwave background and the large-scale distribution of galaxies rule out defects as the sole structure forma- 
■ - - ' tion mechanism, although it is conceivable that some small fraction of the microwave background fluctuations arise 
from defects. In this case as well, direct detection of the gravitational radiation from defects appears improbable. 

The most promising source of a detectable cosmological background of stochastic gravitational waves is a phase 
transition in the early Universe |^,^ . A first-order phase transition proceeds via the random nucleation of bubbles 
of the new phase, which subsequently expand and merge, converting the old phase to the new phase. The coherent 
motion of the bubble walls, which contain a significant fraction of the free energy associated with the phase transition, 
can produce copious gravitational radiation ||^-|^ . The radiation spectrum generically peaks at a comoving wavelength 
corresponding to the Hubble length at the time of the phase transition times the bubble wall velocity in units of the 
speed of light. Remarkably, the horizon scale at the electroweak phase transition falls into the frequency band of 
the proposed Laser Interferometer Space Antenna (LISA) space-based laser interferometric gravitational radiation 
detector and a reasonably strong electroweak phase transition (although much stronger than in the standard 
model) would be detectable with currently planned gravitational wave experiments |^ . 

Besides the bubble wall motions in a phase transition, a related source of gravitational radiation is the subsequent 
turbulent motion of the plasma following the phase transition. Dimensional analysis suggests that turbulence might 
contribute a gravitational radiation background comparable to or larger than that from bubble wall motions [p|,p^. 
In the absence of bubble shape instabilities, the bubbles of the low-temperature phase will expand spherically until 
encountering other expanding bubbles. After the bubbles collide, a region of complex, turbulent plasma motions 
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will result since large amounts of energy are being injected on a particular characteristic length scale. As the phase 
transition completes, the bubble wall motions sourcing the turbulence cease to be effective, and the turbulence damps 
away with a characteristic damping time scale depending on the plasma viscosity. If the bubbles are unstable to 
distortions of their shapes, then the expansion of the non-spherical bubbles can also create additional turbulence. If 
the turbulence is strong, with velocities some non-negligible fraction of the speed of light, significant gravitational 
radiation can be generated during the interval between the initial bubble collisions and the damping of the turbulence 
after the completion of the phase transition. 

In this paper, we quantify these claims by computing the gravitational radiation resulting from an idealized turbulent 
source. We assmne that a source of turbulence exists for some specified length of time, injecting energy on a particular 
length scale at a particular redshift. We model the resulting turbulence as having a Kolmogoroff energy spectrum. 
Details of the turbulence model and discussion of the validity of various assumptions are presented in Sec. II. We 
then compute the generated gravitational waves using the turbulent plasma motions as a source to the wave equation 
(Sec. Ill); the results are then converted to present-day amplitudes and energy densities as functions of frequency. 
Section IV derives the additional gravitational radiation generated by turbulence-induced magnetic fields, showing 
that the peak amplitude from this source will be far smaller than the peak amplitude from the turbulence itself, though 
at a higher frequency. In Sec. V, we apply the results to a generic model of first-order phase transitions, including 
a brief review of hydrodynamic bubble evolution. Section VI discusses the detectability of the resulting backgrounds 
with planned and envisioned experiments. Throughout the paper we employ natural units with c = ?i = fcs = 1. 

A substantial literature on cosmological turbulence appeared three decades ago, when turbulent vorticity was 
considered as a mechanism for initiating galaxy formation JO] . While this particular idea soon fell out of favor due 
to inconsistency with the microwave background isotropy [12| and nucleosynthesis [l3| , some formal aspects of these 
treatments are relevant for this work; see, e.g., p^~p^, which develop phenomenological descriptions of cosmological 
turbulence similar in spirit to that presented in this paper. The hydrodynamic equations in an expanding Universe 
were derived through a transformation of the nonexpanding case in , a special case of a more general theorem [ p^ . 

We emphasize that our results are independent of the nature of the turbulence source. While first-order phase 
transitions are the only obvious source of strong turbulence in the early Universe, the calculations presented here are 
equally applicable to any other potential source of turbulence (see, e.g., ]l9t). 



II. MODEL ISOTROPIC TURBULENCE 

The theory of turbulence was originally formulated over sixty years ago [po| , pT| . But the complexity of turbulent 
motion makes any analysis beyond basic scaling considerations and dimensional analysis intractable. Model isotropic 
turbulence is experimentally tested via wind tunnel measurements on scales small compared to the size of the tunnel, 
and the concepts of a cascade of kinetic energy from large to small scales and the role of viscosity are well established. 
But classical turbulence analysis is done for non-relativistic fluid velocities and incompressible fluids. Here we need 
to model turbulence in a radiation-dominated plasma, potentially with moderately relativistic fluid velocities and 
complications like shock formation. While the theory of turbulence in highly relativistic plasmas is not well understood, 
we will simply extend the nonrelativistic results in the naive manner with the understanding that some corrections 
might apply. 

Consider an event in the early Universe, presumably a first-order phase transition, which converts an energy density 
Acpvac hito kinetic energy of the primordial plasma in some characteristic time scale Tgtir on some characteristic source 
length scale Lg. Here pvac is the total free energy density liberated and k is an efficiency factor which accounts for the 
fraction of the available energy which goes directly into kinetic, as opposed to thermal, energy. The length scale Ls 
must be connected to the Hubble length ~ mpi/T^, which is the only cosmological length scale at early times; 
we write Ls = jH~^. Here T* is the temperature of the Universe when the event takes place and mp\ is the Planck 
mass. Under suitable conditions discussed below, a turbulent cascade will develop in which energy will be transferred 
from larger to smaller scales as eddies of progressively smaller sizes are formed from larger ones. The cascade stops 
at a damping scale Ld when the fluid kinematic viscosity v diffuses the turbulent velocities at the same rate as they 
are replenished from larger scales. We assume that for scales L in the range Ld < L < Ls (the inertial range), the 
turbulence is homogeneous and isotropic. We also must know the enthalpy density w = p + p of the (nonturbulent) 
plasma, which appears in the stress-energy tensor. In our simplified model, any turbulent source in the early Universe 
is determined completely by the physical quantities Pvac, Tstir, -^5, T*, w, and i^. These quantities in turn determine 
Lu, the damping scale, and r, the total duration of the turbulence. Note that a given cosmological model determines 
w and I' from the temperature T*. We also define the wave numbers ks = 2tt/Ls and ko = 2t:/Ld corresponding to 
the largest and smallest turbulence scales. 

The turbulent energy in the cascade is characterized by the stationary Kolmogoroff spectrum 
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where pturb is the kinetic energy density of the turbulent motions. The KohnogorofF constant Ck is of order unity and 
e is the energy dissipation rate per unit enthalpy given by |23] 



e = 2v dkk^E{k) (2) 

Jks 

where v is the kinematic viscosity of the plasma. This spectrum holds for a constant rate of energy flow from larger 
to smaller scales; the amplitude is fixed by the rate of energy dissipation. For a non-relativistic plasma, the enthalpy 
density w is just the mass density of the plasma, while for temperatures large compared to the masses of particles 
in the plasma or for any radiation-dominated plasma, w is 4/3 times the thermal energy density of the plasma. 
Combining the above two equations and solving for the energy dissipation rate gives 

e-f^-^ (3) 

assuming Ck — ^ and ks <C kr). However, E{k) is not yet determined since we do not know the wave number ko 
corresponding to the smallest-scale turbulent motions. 

Before completing the specification of E{k) in terms of the physical variables defining the phase transition, consider 
the time scales involved in the turbulence. Assume that the only peculiar velocities present are (relativistic) turbulent 
velocities with spatial distribution u(x); we employ the Fourier convention 

u(k) = ly dxe*-u(x). (4) 

We retain the fiducial volume factor V to insure consistent dimensions for all quantities; all physical results will be 
independent of V . A statistically isotropic and homogeneous velocity field of an incompressible fluid has the two-point 
correlation function 

(u.(k)^*(k')) = ^P.,(k)F(fc)5(k~k'), (5) 

where 

Pijik) = dij - kikj (6) 

is a projector onto the transverse plane: 

PriP.k = Pk, Pfk, = 0- (7) 

The angular brackets in Eq. (|^) mean a statistical average when the velocities are considered as random variables (see 
Ref. Volume 1, for a detailed discussion). If the fluid is compressible, a second arbitrary function appears in the 
correlation function, proportional to fcifcj, describing longitudinal motions. A specific model for isotropic turbulence 
consists of specifying the function P{k); we assume the power spectrum is a power law, P{k) = Ak", where the 
normalization A and the spectral index n can be deduced from the Kolmogoroff spectrum. The mean square velocity 
of the fluid at any point in space is given by 

(u2(x)) - ^ I dk2P{k) = ^ 1^'" dkk^P{k). (8) 

But this quantity is just the kinetic energy density per unit enthalpy density of the fluid; thus we derive the connection 



E{k) = -^k^P{k). (9) 
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For the case of a Kolmogoroff spectrum, Eq. (hi) implies that 



P{k) ~ l7r2£2/3fc-ii/3 (10) 
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We are interested in the characteristic eddy velocity on a given scale L. From the slope of the KolmogorofF spectrum 
and Eq. (||), it follows that the total turbulent velocity at a given point is dominated by the eddy velocity on the 
largest scale. We can thus estimate the characteristic eddy velocity on the scale L by cutting off the integral in Eq. (||) 
at a wave number — 27r/L corresponding to that scale: 
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(27r)^i/3(^^)i/3_ 



We can also estimate an eddy turnover time scale (known as the circulation time) on a length scale L as the ratio 
of L to the physical velocity vl = + u\Y^'^- We argue below that the physical velocity will be approximately 

bounded by the sound speed of the fluid; for a radiation-dominated plasma, this condition is vl < l/\/3- Making 
the simple approximation that vl = ul until the sound speed is reached, after which time vl is the sound speed, the 
circulation time is 

[ Lv3, otherwise. 

Now the remaining undetermined quantity in the turbulence spectrum, ko, can be fixed via energy considerations. 
Two different cases must be considered separately, depending on whether the duration of the turbulent source Tgtir is 
long or short compared to the eddy turnover time scale ts on the characteristic length scale of the source Ls- First 
consider the simpler case where Tgtir ^ ts- Fully developed turbulence is established in a time on the order of ts, so 
this case gives approximately a stationary source lasting for a time r = Tstir- To keep the turbulence stationary, the 
energy dissipation rate must equal the mean input power of the source: 

e=1P^. (13) 

This expression immediately determines the amplitude of the Kolmogoroff spectrum, Eq. ([^), and comparing with 
Eq. (I) gives 
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Thus the turbulent gravitational wave source is completely determined for this case. The circulation time scale on 
the scale of the source is approximated by combining Eqs. (^2|) and (|l^) to give 
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SO the condition for this case to be valid becomes 
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r.tir » Ls . (16) 



Finally, the Reynolds number for this turbulence is given by 



' ; \ 4/3 r, / 1 \ 4/3 / r4 \ 1/3 

^11) ^) ^ m 



ks J 3 \2tT J \ I^^TstirW 

The critical Reynolds number for the onset of stationary turbulence is around 2000. Early Universe phase transitions 
will generally have Reynolds numbers exceeding this value. 

The alternate case, for Tstir ts-, is more subtle. Here, an impulsive force is imparted to the plasma, resulting in 
a total kinetic energy density equal to the total free energy density of the phase transition times the efficiency factor 
K, coherent on the length scale Ls- The efficiency factor depends on the mechanical details of the stirring process 
and will be a function of mean input power Pvac/Tstir- A cascade of kinetic energy to smaller scales will occur, but 
stationary, isotropic turbulence will never develop because the plasma is not continually being stirred by the source. 
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We can estimate the time for which significant kinetic energy on a given scale lasts. On the largest scale Lg, the 
kinetic energy will last for a time set by the dissipation time scale, approximately equal to the eddy turnover time 
T5. As in fully developed turbulence, this kinetic energy will cascade to smaller scales. The eddies on the largest 
scale will act as a source for eddies on a slightly smaller scale L for a time T5. On the smaller scale, we assume the 
plasma has no kinetic energy at the moment of the impulsive force but rather acquires kinetic energy only from the 
cascade. The smaller-scale eddies are spun up in a time corresponding to the circulation time on the smaller scale ; 
these eddies will last until the large-scale source becomes ineffectual and then will dissipate also on the circulation 
time scale t^. So by this argument, the eddies on a smaller scale L will exist for the same total amount of time as 
the eddies on the largest scale L5, although their establishment and dissipation will be displaced to a slightly later 
time compared with the largest-scale eddies. The same reasoning can then be applied to eddies at successively smaller 
scales, with the following conclusion: on any given scale between Lg and Lj^^ eddies will exist for a total time Tg. 
The only assumption required for this conclusion is that the time scale for establishing eddies on a given scale via 
the cascade from larger scales is the same as the time scale for dissipating the same eddies via the cascade to smaller 
scales. 

The time displacements of the time intervals for the existence of eddies on different scales are essentially irrelevant for 
the generation of gravitational radiation, leading only to some relative phase shift between the gravitational radiation 
at two different frequencies. Therefore, for the purposes of modelling a gravitational wave source, we assume the 
plasma motion consists of kinetic energy simultaneously on all scales within the inertial range, lasting for a total time 
T = Ts, the circulation time on the scale of the turbulence source, with a kinetic energy density spectrum given by 
the Kolmogoroff spectrum, Eq. To normalize the spectrum, we simply treat the total free energy density as being 
injected continually over the time ts rather than as an impulse. Now this case looks just like the previous one, except 
that Tstir must be replaced by rg in Eqs. (O) and (p^): 
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Combining the two cases gives the simple expressions 

\ 27v^TW J 2 \ Kpvac / 

valid for either case, where 

T = max(rs, T^th). (20) 

The eventual expression for the gravitational wave amplitude is only very weakly dependent upon r, so the distinction 
between the two cases is largely unimportant for our results. 

When computing the gravitational wave signal, we will encounter unequal time velocity correlators of the form 

{u,{\<i,t)u*{]^ ,t')) = i^P,^.(k)^(fc,i-t')'5(k-k') (21) 

[cf. Eq. (^)]. The dependence of the function F only upon the time difference t ~ t' follows from the assumption 
that the turbulence can be treated as stationary, with F(fc,0) = P{k). No general form is known for F. However, 
general physical considerations imply that F must be a decreasing function oi t — t' , and we assume that the decay 
of F should have a characteristic time scale on the order of the circulation time on the scale L = 2-K/k. We actually 
will only need to guarantee that F goes to zero no faster than the light-crossing time of L, which is guaranteed by 
causality. 

We have sidestepped the issue of relativistic versus nonrelativistic turbulence. The Kolmogoroff model of turbulence 
phenomenology has only been formulated and tested for turbulence with nonrelativistic velocities in plasmas with 
nonrelativistic equations of state. No general model exists for the opposite situation of a relativistic plasma with 
relativistic velocities. The plasma in our case will also be compressible, contrary to the basic assumption above. For 
a large enough input of energy, plasma velocities may be driven past the sound speed, leading to shock formation. 
We conservatively assume that the sound velocity represents an upper limit to the turbulent plasma velocity, because 
shocks will result in significant thermal dissipation. Note that to the extent that shock fronts retain kinetic energy, 
our ultimate gravitational wave background will be increased relative to the estimates made here, in the case of highly 
relativistic fluid velocities. 

To summarize, this model for cosmological turbulence requires (i) K/Ovac, the energy density converted to turbulent 
motion where pvac is a characteristic energy density and k is an efficiency factor; (ii) L5, the characteristic length 
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scale of the source producing the turbulence (the "stirring scale"); (iii) Tstir, the duration of the source producing 
the turbulence; (iv) T*, the temperature of the Universe at the onset of the turbulence, which in turn determines w, 
the enthalpy density, and v, the kinematic viscosity of the plasma. The assumption of stationary homogeneous and 
isotropic KolmogorofF turbulence then specifies in terms of these quantities (i) the normalization of the turbulence 
power spectrum, (ii) the length scale Lu at which the turbulence is dissipated by viscosity, and (iii) the circulation 
time for any particular turbulent length scale between Ls and Ljj. 

We have neglected the expansion of the Universe in this description of turbulence. If the duration of the turbulence 
r is longer than the Hubble time H~^, then the expansion will produce additional damping of the turbulence as the 
energy density is redshifted. Furthermore, if the circulation time on the stirring scale rg is comparable to or longer 
than the Hubble time, the expansion damping may inhibit the establishment of a turbulent cascade. Particular cases 
should be checked individually, but in general, if a phase transition is strong enough to drive turbulence producing 
an interestingly large gravitational radiation amplitude, it will last for a time short compared to the Hubble time and 
expansion damping will be negligible. This claim can be quantified using the expressions derived in Sec. M below. 



III. GRAVITATIONAL RADIATION FROM TURBULENT PLASMA 



A. General Considerations 



The source of gravitational radiation is the transverse and traceless piece of the stress-energy tensor of a given 
system. For turbulent plasma, the relevant stress-energy tensor is given by 

Tij{x.) ^ wui{x)uj{x.). (22) 

The above expression drops the diagonal (trace) component of the stress-energy because it cannot source any gravita- 
tional radiation. To simplify the problem, we assume (conservatively) that the enthalpy density w remains constant 
throughout space, while the variation of the velocity vector describes the turbulent motions of the plasma. If this 
assumption does not hold, the resulting gravitational wave amplitude will increase. In Fourier space, the stress-energy 
is then given by the convolution 

T,,{k,t) ~ J dqu,{q,t)uj{k~q,t). (23) 

Gravitational radiation is produced by the transverse and traceless piece of the stress-energy tensor. Given an 
arbitrary stress-energy tensor in Fourier space, Tij(k, t), the portion sourcing gravitational radiation can be obtained 
by applying a projection tensor (see, e.g., [p3[): 



Once the source is specified, the gravitational wave metric perturbations hij obey the wave equation 

fh^^2dadh^^ P ^ STrGa^H,, (25) 
a?7^ a arj ai] 

where 77 is conformal time, k is the comoving wave number, and a is the scale factor of the Universe. Note that we 
have defined the tensor metric perturbation as dgij = 2hij. For relevant phase transitions, the duration of the source 
will be short compared to the Hubble time, which means the expansion of the Universe can be neglected during the 
generation of the waves. We can thus drop the expansion drag term in Eq. (|2^) and change variables to physical time 
and physical wave number, obtaining the simple oscillator equation 

hij{k,t) + fc2/iy(k,i) = 87rGny(k,t), (26) 

where dots denote derivatives with respect to t. From this point on, all wave numbers will refer to physical, not 
comoving, quantities. 

The source considered here turns on at a specific time and we assume no gravitational radiation exists prior 
to this time. The initial conditions for Eq. ( p6| ) are simply /iy(k, i*) = /iy (k, i*) — 0. In the Euclidean space 
approximation we have made, the radiation generated cannot depend on the particular value of ts,, so for convenience 
we set <* = in this section. Of course, once the results are translated back into expanding spacetime, the time of 
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the phase transition fixes the energy and length scale associated with the phase transition. The Green function for 
the homogeneous equation is simply 



fcSm 



with G = G = Oatt = 0. The general solution for the wave amplitude is then 

h,j (k, t) = ^ f e{t- t') sin [k{t - t')] n„ (k, t')dt' (28) 
Jo 

where is a step function. 

B. Time Averaging Technique 

Since turbulence is a stochastic process, we cannot compute the exact gravitational waveforms. Our goal is to 
compute the average power spectrum or characteristic amplitude of the waves. We are concerned here only with the 
power spectrum, so consider the quantity 

{h,,ik,t)h*^ik',t)) =i^!^j(k-k') 

x^^^ dh dt2eit-ti)e{t-t2)sm[k{t-ti)]sin[k{t-t2)]n,j{k,h)n*j{k,t2)Y (29) 

The delta- function factor is guaranteed by statistical isotropy of the gravitational waves; we have written this depen- 
dence out explicitly and then changed all factors of k' to k within the angular brackets. To make further progress, we 
need a practical way to deal with the averaging process. We are assuming a stationary, homogeneous and isotropic 
source, so we make the simple assumption that the statistical average can be estimated by either a time or space 
average. To evaluate Eq. (j2^) we use a time average, since all of the time dependence is in the Green functions and 
not in the source terms. Then we have 

{h,,(k,t)h*^(k',t)) =5(k-k')^^^ J dh dt2n,,ik,h)n*^ik,t2) 

dte{t-ti)e{t-t2)sm[k{t~ti)]sm[k{t-t2)], (30) 



X — 
T 

where s is some arbitrary time when the source is active, and T is an interval of time long enough for the average to 
be approximated by the time average. In practice, this will be some time on the order of a few circulation times on a 
given scale. As ti or t2 approaches r, it will not be possible to choose T large enough for a rigorously valid average, 
but this will not appreciably affect our estimates since we are considering only statistical averages for the source 
terms: the time integration is a convenient device for approximating the effect of this averaging, and the averaging 
itself becomes only a rough approximation for durations shorter than the circulation time on a given scale. Since 
we are assuming a stationary source, Eq. (|3^) must be independent of the chosen value of s. We choose an s which 
eliminates the step functions from the integral, keeping in mind the above discussion. 
The integral over t is now elementary: 

1 1 1 

— / dt sin [k{t - ti)] sin [k{t - t,)] - cos [k{t2 - ti)] sin(Tfc) cos [k{2s + T-ti- is)] ■ (31) 

T Js 2 2Tk 

We neglect the second term with respect to the first since Tk ^1: k~^ will be on the order of the light crossing time 
for a given scale, while T will be at least as long as the circulation time on the given scale k~^, so the comparison 
will be valid on all scales except for possibly the largest, where at least the simple inequality Tk > 1 will hold. Since 
the terms are both oscillatory, the comparison really only applies to the size of the prefactors, but this is sufficient 



for our purpose. Now substituting Eq. (31) into Eq. ( pOD and making the substitution y ^ t2 — ti gives 

{h,,{k,T)h*^{k',T)) ^S{k-k')^^^ J^^ dy cos{ky) J\hll,,ik,h)n*jik,h+y). (32) 
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We now use the ti integral as an estimator for the statistical average of the sources, giving 

(/i.,(k,T);i*^.(k',T)) c ^^^^J^ £ dy cos(fcy) (n,, (k, t,)U*^ {k',t, + y)) , (33) 

where the delta-function has been reabsorbed into the statistical average. Note that the average on the right side is 
independent of ti since the source is assumed to be stationary. 

We now have an expression involving the average source correlation at different times, integrated against an os- 
cillating function. Note that the total value is proportional to r, the duration of the source, as it should be for an 
incoherent source. To make further progress, we require a more explicit form for the source average. 



C. Evaluation of the Source Average 



We have expressions for averages of the fluid velocities in the turbulent source; we need to connect these with the 
particular average required in Eq. (33). Writing out the projectors in Eq. (E4) gives 



(%(k,<)n*^.(k',t + y)) = [P,4k)P,5(k) - ^P.,ik)Pabik)] [P,,(k')P,d(k') - ip,(k')Fcrf(k')] 

x{Tab{k,t)T:^{k',t + y)). (34) 

We need to evaluate the expectation value of the stress tensor product. Equation ( p3| ) shows that this product will 
involve the expectation value of four velocity vectors evaluated at two different times. No general solution is known for 
such expectation values for turbulent flow. The simplest (and most conservative) assumption is that the correlation 
function factors into products of pairs of velocities, as for a Gaussian field. Then Wick's theorem applies and we have 



{Tab{k,t)T:^ik',t + y)) 



(27r)6 



w / dqds 



(ua(q, t)ul{q - k, t)) (wc(-s, t + y)u*^{k! -s,t + y)) 
+ {ua{(l,t)ul{s,t + y)) {ubik- q,t)u*a{k' -s,t + y)) + {uaiq, t)ul{k' - s,t + y)) (uf,(k - q, tX(s, i -|- y)) . (35) 

This expression can be simplified using the correlation functions in Eqs. (^) and (^T]), giving 

{Tabik,t)T:^ik',t + y)) ^w^S{k-k') J dq[Paciq)Pbd{^^n) + Pad{ei)Pbc(M^ci)] F{q,y)F{\k-q\,y). (36) 

The first of the three terms in Eq. ( p5| ) does not contribute, since it is nonzero only for the constant offset mode with 
k = k' = 0. After substituting the explicit form for the projectors, Eq. (|6|), setting k = k' from the delta function, 
and simplifying the contractions, we obtain 



(Hy (k, t)n*^ (k', t + y))^ w'S{k - k') J dqF{q, y)F{\k - qU)(l + + P^), 

where we have defined the auxiliary quantities 7 = k • q and = k • k — q. 

Substituting this simple form for the unequal time source correlation into Eq. (|33|) gives 



(37) 



(%(k,T)/j*^.(k',r)) = 



2fc2 



-^(k-k') J dqil + j')il+p') dycos{ky)F{q,y)F{\k-qly). 



Now F{k, 0) = P{k) so we make the further assumption that F can be separated as 

F{k,y) = P{k)D{yk^/'); 



(38) 



(39) 



that is, we have assumed a universal form for the time decay for all k values, with the time argument of F scaling 
with the circulation time on the length scale 27r/fc, and D is some monotonically decreasing function of its argument. 
This is likely a reasonable assumption for fully developed turbulence. On the other hand, we are only concerned with 
the time dependence to the extent that it is integrated against the oscillatory function cos(fct) in Eq. (|3^). Since F 
or D is everywhere positive, the integral itself is oscillatory. If our crude turbulent model were exact, the induced 
power spectrum of gravitational waves would exhibit oscillations. But this is an artifact of the assumption that the 
turbulence begins and ends at precisely defined times. For the present task of estimating characteristic amplitudes for 
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a realistic turbulence source, we instead approximate the time integral by its root-mean-square value. The cos(fcy) 
term will always oscillate on a time scale shorter than the characteristic time for D{yk^^^), as seen from a simple 
comparison of the circulation time to the light crossing time for a given scale L. Thus regardless of the particular 
time dependence of D, we approximate 



^ dycos(fcy)F(g,y)F(|k-q|,y)~^ cos(fcy)P(<7)P(|k - q|) ~ ||p(q)P(|k - q 



(40) 



This approximation replaces the time-dependent function D by the constant D{0). Actually D will decrease with 
time. This will increase the mean value of the integral unless the characteristic time scale for the decrease of D is less 
than k~^, which we have argued will never be obtained, so the approximation in Eq. ( ^ is actually a conservative 
one. 

Substituting this result into Eq. (^sj) and replacing 7^ and by their average values of 1/2 over the integral gives 
the simple approximate form 

{h,,{Kr)h*^{k\r)) ^'^^\fj^'"' s{k--k') J dqP(q)P(|k-q|). (41) 

We now have an expression which can be evaluated for the particular turbulent power spectrum to give the final 
expression for the power in gravitational radiation in terms of the turbulence parameters. 



D. The Power Spectrum 

For a power law power spectrum, the remaining integral in Eq. (^) is elementary. Using the general form P{k) = 



I 



dqP((?)P(|k- q|) = 27rylM dqq 



= 47ryl^ 



d7 (fc^ + q^ - 2kq-f) 

i,2n+3 unun+d 



n/2 



(n + 3)(2n + 3) 2n + 3 n + 3 



(42) 



For the specific case of Kolmogoroff turbulence, the power law is rt = —11/3; then the last term in Eq. (^) is 
dominant. Keeping only this term and inserting Eq. (nm for the power spectrum gives 



(43) 



To make contact with measurable quantities, we evaluate the real-space correlation function 

(/i.,(x,T)/l,,(x,T)) =^ J dkdk' e^(^'-^y- {h,,{KT)h*^ik',T)) 



27V27r' 



G2™2e4/3fc-2/3 / ^fcfc-14/3_ 

2 ' Jks 



(44) 



Now we need to convert this expression to one involving the gravitational wave frequency /. The frequency is 
determined by the scale of time variation corresponding to the spatial Fourier mode k, the circulation time tl given 
in Eq. (^2|). Writing f = and changing variables in the k integral gives 



fs 



(45) 



the numerical prefactor is about 0.007. We define the characteristic gravitational wave amplitude hc{f) per unit 
logarithmic frequency interval (following Maggiore p3]) via 



(/l^j(x,T)/l,y(x,T)) 



#.2 



Kif). 



(46) 
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Note that Eq. (|4^) is smaller than the corresponding expression in Ref. |2J] by a factor of 4 since the tensor metric 
perturbation in Ref. is defined as 5gij = hij whereas ours is 6gij = 2hij [see comments after Eq. (^)]. Comparing 
with Eq. (|5[) gives 

K{f) = 0.12G£-^/W/2/5 (47) 
for frequencies between fs and /d, with the frequency at the stirring scale 

fs f'/'L,'^'. (48) 



E. Relic Gravitational Radiation 

The above expressions apply to the waves generated at the time of the phase transition. We then stretch the 
waves with the expansion of the Universe: the frequency and amplitude are both inversely proportional to the scale 
factor. The latter follows from the fact that the total energy density in gravitational radiation scales like with the 
expansion, and the energy density is proportional to (hh). For turbulence generated at a time when the temperature 
of the Universe was T* , the ratio of the scale factor then to the scale factor now is 

^ = 8.0 X 10- f I^V'VH^) (49) 



ttQ \ 9* J \ T^, 

where is the number of relativistic degrees of freedom at the temperature T, . The Hubble parameter at this time 
is 

SttG ^ STT^g^^ 
with mpi the Planck mass. This gives the relation 



1.65X IQ-^Hz I ^ ) ( )(^) (51) 

\hJ VlOOGeVy VlOO/ ^ ' 

where /* is a radiation frequency at the cosmic temperature T* and / is the corresponding frequency of the radiation 
today. Scaling Eqs. ( ^7| ) and ( ^8| ) by the expansion of the Universe and substituting w = 4pi.ad/3 and e — Kpvac/ (wt) 
gives 

2/3 / r Y^^^ f Ls /lOOGeV\ /lOOV^^ / 



for the characteristic amplitude, which holds for f > fs, and 

-1/3 / 7- \ -2/3 



/s = l.lxlO-Hzf^)"7^ J^] (— ) ■ (53) 

Equations ( ^2|) and (^3|) are our fundamental results. Converting to the characteristic energy density in gravitational 
radiation via the relation 



Kif) = 1.3 X 10-i« (^jj V^r!Gw(/>2, (54) 

where h is the current Hubble parameter in units of 100 km/s Mpc~^ and flcwif) is the energy density in gravitational 
waves per logarithmic frequency interval in units of the current critical density, gives 

. 2. . :o-' {^f (i^)' (a-'" {iY'^- (») 
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IV. GRAVITATIONAL RADIATION FROM INDUCED MAGNETIC FIELDS 



In addition to the turbulent motions, gravitational radiation also may be generated by magnetic fields arising from 
a turbulent dynamo mechanism: generically, the turbulence will exponentially amplify any seed magnetic fields until 
the field strength saturates at equipartition with the turbulent kinetic energy. The characteristic e-folding time scale 
on a give length scale L will be simply the circulation time tl- The mechanism of seed field generation is not clear, 
but seed fields might naturally arise during a phase transition due to bubble wall instabilities combined with surface 
charge densities on the bubble walls and magnetohydrodynamic amplification psf . Once a magnetic field is generated, 
the high conductivity of the primordial plasma will keep the field frozen in. 

It is reasonable to suspect that such a field may give a significant background of gravitational radiation: since the 
magnetic field has a nonzero stress, it will provide a coherent source term in Eq. (|25|). Such a magnetic field will act 
as a gravitational radiation source from the time of the phase transition until the field is damped (or until matter- 
radiation equality, if the field lasts that long), rather than just during the brief period of turbulence. The following 
calculation, however, shows that induced magnetic fields produce a maximum characteristic amplitude of gravitational 
radiation which is always much smaller than the maximum amplitude from the turbulence which generated them. 
The magnetic field gravitational radiation peaks at a much higher frequency, though, and can have a larger amplitude 
than the turbulence-induced gravitational radiation at that frequency. As in the previous section, quantities below 
are physical, except for comoving quantities denoted with a tilde. 



A. General Magnetic Field Considerations 



First, we assume the turbulence-induced magnetic fields are generated almost instantaneously during the time 
of the phase transition. To a good approximation, the turbulence-induced magnetic fields are generated within a 
tiny fraction of the Hubble time , thus we can normalize the magnetic field power spectrum at the time of the 
phase transition. Second, we assume the turbulence-induced magnetic fields are saturated at an equipartition value 
up to a physical scale Lb at the time of the phase transition. We will leave the ratio between the magnetic field 
physical saturation scale to the turbulence stirring scale, i.e. Lb/Ls^ as a free parameter in our final expressions. 
This will make the comparison with the previous turbulence results easier. Generally, we expect Lb/Ls to be on 
the order of 0.003: the turbulence circulation time scales like L^/^, so (^Lb / LsY^'^ gives the ratio of e-foldings of the 
magnetic field on the scales Ls and Lb- Conservatively estimating that the turbulence lasts for a single circulation 
time on the largest scale, a range of Lb/Ls between 0.003 and 0.0017 gives a range of magnetic field amplification 
factors between 10^" and 10'^'^. The exponential amplification makes this estimate robust: making the seed fields 
smaller by a factor of 10^" only reduces by a modest fraction. Third, we assume the turbulence-induced magnetic 
fields are just frozen into the plasma and retain the form of the spectrum until they are damped away by neutrino 
viscosity. Damping of magneto-hydrodynamic (MHD) modes by neutrino viscosity is most efficient before and around 
nucleosynthesis (T ~ 0.1 MeV). At the time of neutrino decoupling (T ~ 1 MeV), the neutrino physical mean free path 
{lydcc ~ 10"'^^ cm) and the Hubble length [H~^^^ « 5 x lO"'^" cm) are comparable, hence all the subhorizon magnetic 
perturbations generated during the electroweak phase transition will be damped away by the time of nucleosynthesis 
(see, e.g., ||l^,^). We do not consider any kind of inverse-cascade mechanism that will transfer small-scale magnetic 
fields to larger scales. Invoking an inverse cascade will spread the magnetic energy to scales larger than Lb and reduce 
the overall gravity wave amplitude. This will also push the gravitational radiation frequencies to smaller values than 
those obtained below. 

A statistically homogeneous and isotropic stochastic magnetic field has a two-point correlation function given by 
Eq. (^ with a power spectrum we denote PB{k). We assume that the turbulence-induced magnetic field exists on 
scales between the saturation scale Lb and the turbulence damping scale Lu. The mean-square value of the magnetic 
field is [see, e.g., Eq. (2.7) of Ref. 

V 



= - dkk^PB{k). (56) 

^ JkB 

We now normalize PB{k) using the fact that the turbulence-induced magnetic field energy density is half of the 
turbulent kinetic energy density 

1 

z on 



Using Eqs. (H) and (|5q), we obtain 
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4^3^;g2/3fc-ll/3 

V 



(58) 



using Eq. (|T^). 

In Fourier space, the turbulence-induced magnetic stress-energy tensor is given by the convolution of the magnetic 
field [see Eq. (2.9) of Ref. l27 



V 1 
(27r)3 47r 



dq 



1 



(59) 



where the explicit dependence is to remind ourselves of the assumption that the turbulence-induced magnetic 
fields are generated almost instantaneously during the time of the phase transition. In addition, we have neglected 
the induced electric field due to the fact that the early Universe is highly conductive. The source for gravitational 
radiation is given by the transverse-traceless projection of this stress tensor, Eq. (p^). 

In the absence of any inverse-cascade mechanism, magnetic fields are just frozen into the plasma and evolve by simply 
redshifting with the Universe's expansion until they are damped away by neutrino viscosity. Therefore, magnetic fields 
act on a longer time scale than the turbulent fiuid velocities. To facilitate the computation, we introduce a comoving 
quantity 11^^- ' (k) corresponding to 11^^^ ' (k, ) via 



n!f^(k) 



(60) 



where k is the comoving wave vector corresponding to the physical wave vector k at the time of the phase transition. 



B. Gravitational Radiation Power Spectrum 

During the radiation-dominated epoch, a oc rj and the homogeneous solutions to Eq. (|2^) are the zero-order 
spherical Bessel functions, i.e. joikrj) and yo{kr]). Defining x = krj and cc* = fcr;,, where ry, is the conformal time 
corresponding to the turbulent source generating the magnetic field, the usual Green function technique yields the 
following inhomogeneous solution for the radiation-dominated epoch: 

Jx, a^W{x') 
where W is the Wronskian of the homogeneous solutions 

W{x) = jQ{x)-^yoix) - yQ{x)-^jo{x) = ^. (62) 

Note that in the turbulence case, the time dependence of the turbulent source is known only statistically. The 
magnetic field, however, is a coherent source, and it evolves by frozen flux until being damped away by neutrino 
viscosity. Therefore in writing down the gravitional wave equation inhomogeneous solution in Eq. (|6^), the explicit 
time dependence of the magnetic source is known and we can immediately perform the time integral, unlike the 
turbulence case. Substituting Eq. ( |62| ) into Eq. (^ij), using the explicit expressions for the zero-order spherical Bessel 
functions, i.e. jo{x) = sinx/x and yo = —cosxjx, and the approximation for the scale factor in the radiation- 
dominated epoch a{ri) ~ iJo?7\/f^rad, we obtain 

~ 87rGn,^fVk) - 

fe|f (k,7?)c. ' ak,V*,v), V<Tk, (63) 

where -qj, corresponds to the conformal time at which the magnetic perturbation comoving wave number k is damped 
away by neutrino viscosity. Here we have abbreviated 

^{k,V*,v) ^ ^{kv*,kTj) = [ c^j^' !!!1[M!L_JL]1. (64) 

Jr,. V 

It is simple to see that ^ is an oscillating function with a monotonically decreasing amplitude of oscillation; the 
amplitude decays more slowly than the 77"^ dependence of free gravitational waves, since the wave is continually 
sourced by the magnetic field. 
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As in the turbulence case, we are interested in the average power spectrum of the waves, so we consider the quantity 

2 



From Eq. (|60|), we have 



SttG 



kriH^n 



rad 



(n|f)(k)n(f*(k')). 



(65) 



(66) 



As in Eq. (2.20) of Ref. the two-point correlation function (n|^''(k, i*)n[^''*(k', i»)) at the time of the phase 
transition can be written as: 



V 

where the tensor structure Aiijim is [Eq. (2.21) of Ref. p^] 



|n(^)(fc,t,)|25(k-k'), 



(67) 



M^jiM = Pu{k)Pjmi'k) + P,;,„(k)P,7(k) - P„(k)Pi™(k) 



(68) 



and satisfies Mijij = 4 and Maim = Alyn = 0. Then using Eqs. (|59|), (|24|), and (^, a similar calculation as in the 
previous section (see also the Appendix of Ref. |E^) gives 



(nlf (k,iOn,(f*(k',i,))^-1^5(k-k') I dqFs(g)PB(|k-q|)(l+7')(l+/32), 



(69) 



where as in Eq. (|3^) we have defined 7 = k • q and /3 = k • k — q. In deriving Eq. (|69|) , we have assumed the 
turbulence-induced magnetic field to be Gaussian, as in the case of the turbulent fluid velocities, and hence we can 
apply Wick's theorem. Comparing with Eq. (^^, replacing 7^ and by their average values over the integral of 1/2, 
and using Eq. (58) gives 



|n(^)(fc, ^ J dqP(g)P(|k - q|). 



(70) 



This integral has already been done in Eq. (|4^), except that now the lower limit for the physical wave number is 
instead of ks; hence 



|n(^)(fc,t,)|2 ^ ^^2^-4/3fc-ll/3fc-2/3^ 



8V 

Equations (|65|), (p6[), (|67|), and ([TI]) together then give 



{h^\k,rj)h\^>*ik',r^)) 



(B)*,r, 



8647r^G2 u;2e4/3fc-17/3^^2/3^28/3 



y2 



^ ^rad 



-e'(fc,^*,^)^(k-k'), 



(71) 



(72) 



where we have approximated a ~ Horjy/ f2rad while the Universe is radiation dominated, and we have converted to 
the comoving quantities V — V/a^, k = ka^,, S{k — k') = S{k — k')/a^. As in the previous section, we evaluate the 
real-space correlation function to make contact with measurable quantities: 

(/i(f)(x,r;)/i|f)(i,r7)) - dkdk' e^^^'-^yHh[f\kv)h';f^*i^',v)) 



i2n) 

M,r2r2,,,2=.4/3r:-2/3^28/3 ^kn 

54nGwe k^ a, [ dkk-''/'eikv*,v) 



^ ^rad 



(73) 
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C. Relic Gravitational Radiation 



As in Eq. (^), we define the characteristic gravitational wave amplitude hi^\f) per unit logarithmic comoving 
frequency interval via 

(/.(f)(i,,)/.|f (X,,,)) ^ ^hi^^\Lv), (74) 

and the statistical average on the left side implies that the time dependence on the right side is not the exact time 
dependence of the gravitational wave but only the time dependence of its amplitude (i.e., the oscillations are averaged 
over). Since Eq. ( |63[) gives the exact time dependence of the gravitational radiation (as opposed to the turbulence 
case, when we only know the statistically averaged time dependence), the standard wave dispersion relation holds: 
k = 27r/. (In contrast, for the turbulence source, we only know the statistically averaged time dependence, so we use 
only the approximate dispersion relation / = .) Then comparing with Eq. (^) gives 

/ o3/2 \ ^,,,^2/3 f-1/3 14/3 

which depends on the function ^ which we define as the amplitude of the oscillations in ^(?y), times the numerical 
factor %/2/2 to convert to a root-mean-square value, in accordance with the definition of the characteristic amplitude 
Eq. (0). In deriving Eq. (^5|), we have used the fact that after the conformal time ?7ond when the magnetic fields are 
damped away via viscosity and cease to be an efficient source of gravitational radiation, the characteristic amplitude 
he {1,11) '^il^ simply scale inversely with a. Note that he (/,'7o) is only weakly dependent on T^end, which occurs 
through the upper limit of the integral in ^(/, ?7*, ?7Gnd). The time r/ond technically depends on the scale considered, 
but for simplicity we simply use the saturation scale Lb at which the gravitational radiation peaks. 

The function ^{kri^,,kr]) cannot be expressed in terms of elementary functions, but it is simple to obtain an upper 
bound. Since the amplitude of the oscillations in is monotonically decreasing, the amplitude at the initial time gives 

ah*j'v)<-^, (76) 

which is useful for constraining the gravitational wave amplitude. 

In an analogous calculation to the previous section, writing Jb = Jbo-* and Jb = Lg^ and approximating _ff* ~ 



i/oV^rada* ^, we obtain 
.(-)(/,.o)..1.9xI0-(^)^^^(^ 



-2/3 , ^5/3 ^iooGeV\ flOOV^' ( f \ ~ 

] ^{f,V*,Vcnd)- (77) 



IB , 



This characteristic amplitude is valid for f > /b, where [using Eq. (]5l|)] 



H-^J VlOOGeV; VlOO/ 



/B = 1.65xl0-Hz 77^ (78) 



The corresponding energy density in gravitational waves per logarithmic frequency interval in units of current critical 
density is 



V. FIRST-ORDER COSMOLOGICAL PHASE TRANSITIONS 



The most likely mechanism for creating turbulence with a large energy density is a first-order phase transition. 
Such a transition is controlled by an effective potential for some quantity which functions as the order parameter of 
the phase transition. Initially, the Universe sits in a minimum of the effective potential. As the Universe expands 
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and cools, the effective potential develops a local minimum at a different value of the order parameter; this new local 
minimum eventually evolves to be the true minimum energy state. Then the order parameter wants to evolve to the 
new minimum. If a potential energy barrier exists between the old local minimum and the new true minimum, the 
phase transition must occur via quantum tunnelling through the barrier or thermal fluctuations over the barrier. As 
a result, bubbles of the low-temperature phase are nucleated at random places in the high-temperature phase. The 
energy difference between the two phases creates an effective outward force on the bubble, causing it to expand. Once 
this outward force from the energy difference balances the inward hydrodynamic force from pushing plasma outwards, 
the bubble reaches an equilibrium and expands at a constant velocity. We will consider only the case of quantum 
tunnelling, applicable to a strong first-order phase transition with a high barrier between the two phases. In this case 
the nucleated bubbles are spherical and negligibly small compared to the horizon scale . The more complex case 
of thermally activated bubbles has been considered in 



29|. 



A. Turbulence 



In general, the rate for nucleating a bubble will be the exponential of some tunnelling action, F oc exp(S'(i)). As 
a simple model of a phase transition, we expand the action S into a power series in time and keep only the constant 
and linear terms. This gives a characteristic bubble nucleation rate per unit volume | pO[ | 

r = Toe^* (80) 

so the quantity sets the characteristic time scale for the phase transition. Numerical calculations show that 
the largest bubbles reach a size of order P~^vi, by the end of the phase transition where Vb is the bubble 
expansion velocity, assuming the bubbles remain spherical as they expand. In general, f3 is expected to be of the order 
41n(mpi/r)i/ ~ lOOH for a Hubble rate H 

A first-order phase transition is generically described by several parameters: (i) a = Pvac/Pthcrmai = 4/5vac/3w, the 
ratio of the vacuum energy associated with the phase transition to the thermal density of the Universe at the time 
(which characterizes the strength of the phase transition); (ii) k, an efficiency factor which gives the fraction of the 
available vacuum energy which goes into the kinetic energy of the expanding bubble walls, as opposed to thermal 
energy; (iii) (3, which sets the characteristic time scale for the phase transition; (iv) Vb, the velocity of the expanding 
bubble walls, which set the characteristic length scale of the phase transition; (v) T*, the temperature at which the 
phase transition occurs. 

Once the bubbles expand and percolate, much of their kinetic energy will be converted to turbulent bulk motions 
of the primordial plasma (for an illustration, see the numerical evolution of two scalar field bubbles in Ref. |^). The 
energy density contained in a bubble wall of radius r scales with r'^, the bubble volume. As the phase transition 
ends, far more small bubbles have been nucleated than large ones, but the energy density in the large ones dominates 
the total energy density We therefore make the approximation that turbulent energy is injected on a stirring 

scale Ls — V},P~^ corresponding to the size of the largest bubbles. The stirring will last for roughly Tgtir — (3~^, the 
duration of the phase transition. The duration t of the turbulence then follows from Eqs. (EG) and (Esh as 



T = (3 ^ max 



3\/2 Vb 
' 2 (Ka)i/2 



(81) 



The fundamental symmetry breaking mechanism which drives the phase transition determines some effective po- 
tential for bubble nucleation. The difference in energy density between the two phases and the bubble nucleation 
rate are both determined by this mechanism. Thus the parameters T*, /3, and a are all determined directly by the 
underlying physics, and are precisely calculable to some given order in the various particle interaction strengths. On 
the other hand, the bubble velocity Vb and the fraction of kinetic energy into the bubbles k depend on the detailed 
microphysics involved in the bubble propagation through the relativistic plasma and are not determined from general 
properties of the effective potential. Generally, the larger the vacuum energy density driving phase transition, the 
higher bubble wall velocities Vb will be obtained. 

The hydrodynamic boundary between a lower-energy phase and a higher-energy one can propagate via two modes, 
detonation and deflagration. Details of these modes in the case of spherical geometry are known . For a detonation 
front, the velocity of the phase boundary exceeds the sound speed in the fluid, so that a shock forms at the burning 
front. In the opposite case, a deflagration propagates slower than the sound speed and piles up an overdensity of 
fluid in front of it, like a snowplow. The boundary conditions for a detonation are more restrictive, so that once the 
energy densities and pressures are specified in each phase, the complete solution for the propagating detonation is 
determined. In this case, we have M% 
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(82) 



and the approximate form |9| 



1 



1 + Aa 



Aa 



4 

27 



3a\ 

"2"; 



1/2 



(83) 



with A = 0.72. If the bubbles propagate as a deflagration front, no such general relations apply. However, it has 
been argued that for relativistic plasmas, instabilities in the bubble shape will accelerate the bubble walls and the 
hydrodynamic expansion mode is unstable to becoming a detonation. For this reason, in the following analysis, we 
will assume Eqs. (^) and (|8^) hold. We also assume that a ^ 1 to simplify further Eqs. (|8^) and (|8^), which will 
generally hold for realistic phase transition models; for unusual cases with very strong detonations and a > 1, the 
following formulas must be corrected. The duration of the turbulence is then given by the second term in Eq. (|8l]), 
becoming 



9/4 



(3- 



-3/4 



(84) 



B. Relic Radiation from the Phase Transition 



The characteristic gravitational wave amplitude from turbulence becomes 

2 / i \ 1/3 



(3 J \ T* J \ g* 




-11/4 



with the characteristic frequency 



fs ~ 5.7 X IQ-^Uza^/* 



1/6 



viooGev J viooy 

The corresponding energy density per logarithmic frequency interval is 

2 



Qcwh^ ^ 2.7 X 10-i"ai5/4 



p) \ioo) I/. 



-7/2 



(85) 



(86) 



(87) 



In a first-order phase transition, the expanding, colliding bubbles are themselves a potent source of gravitational 
radiation For our idealized model phase transition with spherical expanding bubbles, the ratio of the maxi- 
mum amplitude of gravitational radiation due to turbulence to the maximum amplitude due to bubble collisions is 
approximately 



ifs) 



^bub (/"max) 



0.18a-3/8^ 



so only for a < 0.01 will the amplitude of the turbulent signal be larger (although in this case, turbulent damping due 
to the expansion of the Universe is significant and our estimate for the turbulence gravitational wave amplitude may be 
significantly too large). For realistic models with interesting gravitational wave production, the turbulence amplitude 
will be subdominant to the bubble amplitude, but non- negligible. This is in contrast to the naive dimensional estimate 
of the turbulence gravitational radiation in Ref. |^ which gave a somewhat larger value. The frequencies at which 
these maximum amplitudes occur scale differently with the parameters: 



fs 



fn 



l.la3/4. 



(89) 



The different scaling arises because the duration of the phase transition Tstir sets the characteristic frequency for the 
radiation from expanding bubbles, while the circulation time on the stirring scale rg sets the characteristic frequency 
for the radiation from turbulence. 
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Note that the gravitational radiation in the bubble case has a long tail in the amplitude, hc{f) oc f~^^^, while 
turbulence driven at a single scale drops off very quickly like /ic(/) oc /~^^/^. The tail for bubble collisions arises in 
the case of bubble collisions because at any given moment, the characteristic frequency of radiation from the collision 
of two bubbles is wt, /d, where d is the size of the colliding region. Since d ranges from zero to the maximum size of the 
smaller bubble as the bubbles expand, the gravitational radiation is produced over a wide range of frequencies. This 
tail of the frequency spectrum is somewhat model-dependent, and will be modified if the bubbles are not spherical. 
Departures from sphericity could arise from thermal activation over the potential barrier, resulting in non-spherical 
nucleation, or from shape instabilities as the bubble expands. The results for expanding bubbles also depend on the 
thin-wall approximation, namely that the width of the bubble wall is small compared to the radius of the bubble. 
While this approximation will be very good for relativistic detonations, it will not be as good for deflagrations. 

The gravitational wave signal from turbulence from a single stirring scale Lg is somewhat more generic, although 
if the phase transition does not proceed via detonation the specific expressions for k and Vb in Eqs. (|8^) and ( |8^ ) 
will not hold. However, the single-scale assumption obviously will never be exactly correct; any realistic source like 
a phase transition will deposit bulk kinetic energy over a range of scales. The energy density in bubble walls of a 
given size will generically peak at a scale comparable to Vb(3~^ that we have taken for Lg, because the kinetic energy 
in the wall of a bubble of radius r scales like so the energy distribution is heavily weighted towards the largest 
bubbles. Analytic expressions for the size distribution of bubbles, the fraction of space taken up by bubbles, and 
related quantities are given in Ref. On the other hand, the stirring scale appropriate to the collision of two 

bubbles of unequal radius is not entirely clear: some turbulence will clearly be created on the scale of the smaller 
bubble, but since the larger bubble has greater energy density in the wall, a significant part of the energy will remain 
in coherent motion determined by the larger bubble. 

In realistic cases, the gravitational wave amplitude spectrum in Eq. (^) must be convolved over a range of stirring 
scales. A specific model of the distribution of stirring scales in a first-order phase transition is beyond the scope of 
this paper. However, we can make a rough estimate of its effect. Assume that the actual turbulence source stirs 
the plasma over a range of frequencies A/s- The actual bubble size distribution has a significant tail towards larger 
bubbles |^^. If the same total energy goes into gravitational radiation as in the single stirring scale case, then the 
characteristic amplitude hc{fs) will be reduced by a factor of order (/s/A/s)^/^. This very crude estimate neglects 
the strong dependence of the amplitude on the stirring scale and employs only a box-shaped energy density spectrum, 
but the general scaling is correct. Generically, the distribution of bubble sizes in a model phase transition points to 
A/5//5 on the order of a few (see ||30|| ), but a more precise estimate requires a detailed model of stirring in a phase 
transition. As a rule of thumb, when estimating the gravitational radiation background from turbulence arising from 
a phase transition with a single stirring-scale model of the turbulence, the resulting amplitude may be overestimated 
by a modest factor. 



C. Relic Radiation from the Induced Magnetic Fields 



For the magnetic fields from the turbulent dynamo mechanism, the characteristic gravitational wave amplitude 
becomes 

with the characteristic frequency 

/. . 2.9 X 10- Hz (^] f A) f (^) . (91) 



The corresponding energy density per logarithmic frequency interval is 

f]|f^(/>2 ^ 4.8 X (^)"'^' ir ) ^"'(/'?-/^e„d). (92) 



-2/3 



The ratio of the maximum amplitude of gravitational radiation due to turbulence-induced magnetic fields to the 
maximum amplitude due to the turbulent fluid today is approximately 

h^^\fB) 2_7^3/8 ( A.) ( ^V'^"(/B^-/B'7o„d). (93) 
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The ratio of the frequencies at which these maximum ampUtudes occur is 



Ib _ , c 1 „ -3/4 f 



The scahng with a arises because the circulation time on the stirring scale T5 sets the characteristic frequency for the 
radiation from turbulence, whereas for magnetic fields, fs = L^^. 

The value of rjcnd corresponding to the scale Lb can be determined via consideration of the neutrino viscosity 
(see |3^) but ^(Ar;*, fcryond) is only weakly dependent on ?7cnd so we do not compute it here. Instead, we derive an 
upper bound on the amplitude. The approximate relation a*i/* ~ 1/??*, valid during radiation domination, gives 
fBf]* — /Lb- Then Eq. (fz^ ) combined with the frequency dependence in Eq. ( |85| ) gives 

<27a-^V'%J^\'" . (95) 



As discussed above, Ls / Lb — 300 generically, so the peak characteristic amplitude from the magnetic field at frequency 
Jb will always be negligible compared to the peak characteristic amplitude from the turbulence at frequency fs ■ The 
turbulence gravitational waves drop so quickly with frequency, however, that the magnetic field gravitational waves 
will give a larger characteristic amplitude at . 



VI. POTENTIAL DETECTABILITY 



The detectability of a given stochastic background depends on both its characteristic frequency and its amplitude. 
The Laser Interferometer Gravitational-wave Observatory (LIGO) is nearing the commencement of scientific 
observations; it is comprised of two facilities in the United States, each essentially a Michelson interferometer with 
an arm length of 4 kilometers. LIGO has sensitivity to gravitational radiation in the frequency range from 10 to 
1000 Hz. Seismic noise prevents useful gravitational wave detection from the surface of the Earth at frequencies lower 
than about 10 Hz. Cross-correlation of the two LIGO detectors, along with several smaller laser interferometers and 
bar detectors at other sites around the world, allow a clean detection of stochastic signals, since widely separated 
detectors have no correlated sources of noise. Detailed estimates shows that in this frequency range, LIGO will be 
able to detect stochastic gravitational wave backgrounds with a characteristic amplitude of around hc{f) ~ 3 x 10^^^ 
at / ~ 100 Hz after integrating for four months p5|-p^]. These levels will hopefully be obtained within three years. 
Planned technical improvements are projected to reduce this threshold amplitude by another factor of 10 on the time 
scale of a decade. 

The other major gravitational wave observation program, the Laser Interferometer Space Antenna (LISA) ||^, is a 
cornerstone mission of the European Space Agency in partnership with NASA. Current design studies envision three 
spacecraft arrayed in an equilateral triangle with an arm length of around 5 x 10^ kilometers with laser interferometry 
between each of the three pairs of arms; the spacecraft configuration will trail the Earth's orbit by about 20°. LISA 
will likely be sensitive to a frequency range from around 0.0001 Hz to 0.1 Hz. The detection of stochastic backgrounds 
with LISA is more complicated than with LIGO, because any pair of interferometers formed by EISA's arms share one 
arm in common, so it is not possible to cross-correlate two independent interferometers with uncorrelated noise. It was 
originally believed that this limited detection of a stochastic background to the level of the instrument noise power 
because there would be no way to distinguish between instrumental noise and a background signal. This noise level 
corresponds to a stochastic background amplitude of around hc{f) = 10~^^ at 0.01 Hz. It has now been realized that 
if the complete time series data for positions of 6 independent test masses are recorded, so-called Sagnac observables 
can be synthesized which are highly insensitive to various kinds of noise in the system ]39[ |, including one which is 
largely independent of low-frequency stochastic gravitational wave backgrounds, allowing a direct measurement of 
the system noise pO| . This results in a significant improvement in the ability of the system to measure stochastic 
backgrounds ||4l[| . For one year of observation, this kind of analysis could in principle give sensitivities comparable to 
two independent Michelson interferometers, reducing the threshold hc{f) by a factor of [ft^^^ for observation over 
a time t, or hc{f) ~ 4 x lO'^^ at 0.01 Hz over one year of observing. Such a sensitivity level depends on a precise 
understanding of the system noise properties and elimination of other correlated noise sources between the various 
arms of the detector, which is only partially practicable. Flying and cross-correlating two independent LISA-like 
detectors 42,4^] is still clearly preferable for detecting stochastic backgrounds. 

For stochastic background detection at LISA frequencies, raw sensitivity is not the only issue. White dwarf binaries 
in our galaxy will produce an approximately stochastic gravitational wave background which probably becomes 
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comparable to the LISA sensitivity limits for frequencies below about 10~^ Hz Q]. Detection of such a signal will be 
interesting in its own right, but will effectively provide a lower limit of around 10^^ Hz to the stochastic background 
signals which are detectable, until gravitational wave detectors improve to the point of having enough directional 
sensitivity to distinguish sources in the galactic plane from sources distributed isotropically. 

The characteristic gravitational wave frequency for turbulence from known phase transitions is not promising for 
detection in the near future. For the electroweak phase transition at ~ 100 GeV, Eq. ( ^6|) shows that a, the ratio 
of the vacuum energy density to the thermal density at the time of the phase transition, must be of order 0.1 for the 
frequency maximum to be as high as fs = 10~* Hz, if (3/ H^, takes its characteristic value of 100. This frequency 
is the lower limit to what LISA might be able to detect. The amplitude at this frequency for a = 0.1 would be 
hc{fs) — 2.8 X 10"^'^, more than two orders of magnitude smaller than a LISA Sagnac configuration could detect at 
this frequency. Any push towards higher frequencies via a shorter phase transition further reduces the characteristic 
amplitude, since hc{f) oc (i?*//?)^. For an extreme case with a — 1 and /? = lOOOiJ*, the characteristic frequency 
is near LISA's maximum sensitivity, fs — 5.7 x 10~^ Hz, with a characteristic amplitude hc{fs) — 3.8 x 10"^*^. 
This amplitude is an order of magnitude smaller than the LISA Sagnac sensitivity at this frequency. An analysis 
of the electroweak effective potential in a large class of supersymmetric models shows that for models with 

large values of a, generally [3 < 1007?*, and a is never as large as unity Other well- motivated extensions of the 
standard model may result in a very strong electroweak phase transition (e.g., Q). 

Satellite missions to probe stochastic backgrounds to lower frequencies have been discussed , which would involve 
multiple spacecraft arrayed at separations on the order of 1 a.u. Such configurations would be a more natural match 
for the frequency scale of electroweak turbulence, although dealing with the binary foreground signal would still be a 
major hurdle. 

Phase transitions at lower temperatures, like the QCD phase transitions, have larger characteristic length scales 
and thus even lower frequencies for gravitational radiation. Speculative phase transitions could occur at energy scales 
higher than the electroweak scale, resulting in higher characteristic frequencies. However, a higher energy scale also 
translates into a smaller characteristic amplitude, and it is not possible to give a set of parameters with a ^ 1 for which 
cosmological turbulence would be detectable in LIGO. LISA could detect the turbulence from a range of imagined 
phase transitions at energy scales above the weak scale, but at present no compelling theoretical motivation for such 
phase transitions is at hand. 

In contrast to turbulent sources, the expanding bubbles in a first-order phase transition, which drive the turbulence, 
are themselves a strong source of gravitational radiation |^ and are a much more promising source of detectable signals 
from the electroweak phase transition. The difference between the detectability of the two sources is essentially the 
factor of a^^'^ in Eq. (|89|), arising from the different time scales of the sources. The characteristic frequency for the 
expanding bubbles is set by the phase transition time scale (3~^ because the bubbles expand and percolate in this 
time. For turbulence, the time scale is instead the circulation time on the stirring scale. The turbulent fluid velocities 
are significantly smaller than the bubble expansion velocities unless the turbulent flows are near the sound speed, 
giving lower characteristic frequencies. (Extremely strong turbulence with relativistic fluid velocities would likely 
produce gravitational radiation in a more detectable range of frequencies and amplitudes, but the amount of energy 
in turbulent motions is limited by shock formation and heating, and turbulence is not understood in this regime.) 

Our results in this work supercede the dimensional estimates in Ref. |^, which predicted that the turbulence 
signal could be significantly larger than the bubble signal at similar frequencies. Resulting optimistic calculations of 
turbulent signals from the electroweak phase transition, e.g. Ref. |^^, unfortunately do not hold up to more detailed 
analysis. We emphasize, however, that the results presented here apply in a generic way to any turbulence in the early 
Universe, and the search for stochastic gravitational radiation backgrounds in the frequency range from 10~^ Hz to 
1000 Hz is in part a search for unanticipated, dramatic physics at energies above the electroweak scale. Perhaps we 
will be lucky. 
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